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Abstract 

A novel method to derive stationary solutions of the Vlasov-Maxwell system is established. This 
method is based on the assumption that the deviation of the velocity distribution from the Maxwell- 
Boltzmann distribution can be expanded by the Hermite polynomials. By applying our method, a 
new two-dimensional equilibrium is derived, which may provide an initial setup for investigations 
of three-dimensional collisionless reconnection of magnetic fields. 
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I. INTRODUCTION 



It goes without saying that magnetized plasmas play important roles both in nature and 
laboratory experiments. For example, in the solar atmosphere and the terrestrial magneto- 
tail, magnetic reconnection [TJ [2} [3] is considered to be a plausible mechanism converting 
magnetic energy into kinetic energy of particles. Furthermore, in supernovae, gamma-ray 
bursts, and other astrophysical jets, some plasma instabilities develop when the ejecta in- 
teract with the ambient medium [I]. Such plasma instabilities are important also in laser 
experiments. 

There are some approaches to describe the behavior of such plasmas and electromagnetic 
fields. One is the ideal or non-ideal magnetohydrodynamics (MHD). This is based on the 
fluid treatment implicitly assuming the Maxwellian velocity distribution for particles. On 
the other hand, the kinetic treatment of plasmas leads to the Boltzmann equation coupled 
with the Maxwell equations (see, for example, Ref. j5]). By using this approach, we can 
deal with plasmas even if the velocity distribution deviates from the Maxwellian. Neglecting 
the effect of collisions, we obtain the Vlasov-Maxwell system that can describe sufficiently 
rarefied plasmas. 

In particular, the magnetic reconnection is of great interest as an active phenomenon 
of plasmas. This phenomenon, which also occurs in collisionless plasmas of the solar at- 
mosphere and the terrestrial magnetotail, requires an anomalous resistivity, inertia-based 
processes or off-diagonal terms of the electron pressure tensor. Recently, a great deal of 
investigations of the collisionless magnetic reconnection have been performed by theoretical 
analyses U\ U U ED] , particle simulations pH 121 H3 E] , MHD simulations pH EH], and 
a Vlasov simulation |17j . 

Stationary solutions of the Vlasov-Maxwell system provide self-consistent configurations 
of collisionless plasmas and they are the convenient initial setup for investigations of the 
collisionless magnetic reconnection. Therefore, some authors have proposed such solutions 
[IHl [IHl EQl EH [221 [231 EH ESI ETj. The Harris sheet [18] is one of the most well-known 
equilibrium assuming the Maxwellian velocity distribution. This equilibrium is widely used 
as an initial setup for simulations of the magnetic reconnection. The Bennett pinch [19] 
also assumes the Maxwellian velocity distribution and provides an equilibrium in cylindrical 
configuration. The so-called BGK solution [20J treats an electrostatic equilibrium. The 
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effect of the temperature gradient is included in Ref. [2H [221 [27]. Ref. [26] investigates 
two-dimensional case. 

In the present paper, we focus on a method to construct stationary solutions of the 
Vlasov-Maxwell system. By applying our method, we can reproduce some results of previous 
studies. In other words, our method is a generalization of the methods presented in these 
studies: it can treat two-dimensional plasmas with the velocity distribution deviating from 
the Maxwellian. The generalization enables us to explore new solutions. By way of example, 
we construct a new two-dimensional equilibrium, which may provide an initial setup for 
three-dimensional magnetic reconnection studies J2E1 ED] • 

This paper is organized as follows. In SecJTTJ we describe the procedure of our method. 



In Sec, III, some results of previous studies are reproduced. We propose a new equilibrium 



derived by applying our method in SecjTV] We conclude this paper in SecjV} Appendix gives 
some algebraic relations used in this paper. 



II. FORMULATION 



We deal with the Vlasov-Maxwell system. The Vlasov equation describes the kinetic 
evolution of the distribution function fj (t, x, v) of particles j (= i for ions and e for electrons) 
in the phase space (x, y, z, v x , v y , v z ) and the Maxwell equations describe the evolution of 
the electromagnetic fields. Here the cartesian coordinates in the real space are (x, y, z) and 
the corresponding coordinates in the velocity space are (v x , v y ,v z ). This system provides an 
exact model for the behavior of collisionless plasmas (see, for example, Ref. [31]). In this 
section, we derive stationary solutions of the Vlasov-Maxwell system under the assumptions 
described in the following subsection. 



A. Equations 

We consider a stationary plasma. It is assumed to be homogeneous in the z-direction. 
Thus physical variables do not depend on t and z. In this circumstances, the Vlasov equation 
expressed as 

vM- I rj^ ■ '<< (cE x - v z B y ) ^+^~ (cE y + v z B x ) ^ + ^~ (v x B y - v y B x ) ^ = 0, 
ox ay rrijc ov x rrijC ov y rrijC ov z 

(1) 



governs the distribution function fj for particles j with the charge qj and the mass rrij. 
Here c is the speed of light. E x , E y and B x , B y represent the x, y-components of electric and 
magnetic fields, which are functions of x and y. These electromagnetic fields are written as 

d<f> d<p 
Ex = -d~x' Ey = ~dy> (2a) 

R - 9Az B - dA * (2h) 

by introducing the scalar potential <p(x,y) and the ^-component of the vector potential 
A z (x, y). The other components are assumed to vanish. These potentials satisfy the Poisson 
equations; 

W + W = ~ P ' (3a) 
d 2 A z d 2 A z 

~dx^ + ^ = (3b) 

where p(x, y) and j z (x, y) are the charge density and the z-component of the electric current 
density, respectively. These are expressed in terms of fj(x,y,v x ,v y ,v z ) as 



OO POO 



P = / dv x I dv y j dv z fj, (4a) 

- J — oo J — oo J —oo 

POO poo poo 

jz = / dv x / dv y I dv z v z fj, (4b) 



which close the system. Solutions of Equations ([TJ-Q give self-consistent configurations of 
collisionless plasmas. 



B. Derivation of stationary solutions 



In seeking solutions of Equations (JT[)-((4]) , we assume that the distribution function takes 
the following form; 



fj 



y^g jtn (A g )H n [v z j 



Vj ) 



n=0 



n ; 



7f 3/2^3 



exp 



v 2 + v 2 + v 2 



2 Qj 



rrijV 2 



(5) 



where rij and Vj are the average number density and the thermal velocity of particles j, 
respectively. We assume they are constant. H n represents the nth-order Hermite polynomial 
and gj >n is its coefficient, which is assumed to be a function of A z . The distribution function 
in this form describes the deviation from the Maxwellian in terms of the Hermite polynomial 



expansion. The spatial dependence of the distribution function is described through those 
of potentials <fi and A z . 

At first, we substitute the expressions for fields ^ and the distribution function ^ into 
Equation ([!]). Using the recurrence relation (A5), we obtain 

dg j<n 2(n + l)qj 



. n=0 



dA, 



rrijVjC 



~9j,n+l 



H n (v z /vj) 



n 3 

yj-3/2y3 



exp 



2 Qj 



rrijVj 



(6) 

Then we multiply both sides of this equation by H m (v z /vj) and integrate with respect to v z 
to obtain the recurrence formula for <? J)n , 



9j,n+l 



rrijVjC dgj i7l 
2{n + l)qj dA z ' 



(7) 



Here we have used the orthogonality relation of the Hermite polynomials (A2). As a result, 
9j,n{n > 1) are determined once gj^ is given as a function of A z . 



On the other hand, substitution of Equations ([5| and d7j) into Equations Q leads to 

2qi4> \ ( 2q e < 



p = niqig ifl exp 



m,t>„- 



+ n e q e g efl exp 



m e v1 



Jz 



2 dA. 



exp 



2qj(p \ m e n e v 2 e dg efi 
rriivf) 2 dA 2 



exp 



2q e 



m e v1 



(8a) 
(8b) 



where the orthogonality relation (A2) has been used again. Substituting these expressions 
into Equations pj), we obtain 

2qi<f)\ ( 2q e 



d 2 (f) d 2 4 

dx 2 dy 2 

d 2 A z d 2 A z 
+ 



dx 2 



dy 2 



An 
-2tt 



rnqigijo exp 
miUiV 2 dg ifi 



+ n e q e g efi exp 



m e vi 



exp 



2 qi 



m e n e v 2 dg efi 



exp 



2q e 



m e v 2 



(9a) 
.(9b) 



2 dA 7 . \ niivf J 2 dA, 

Thus, once we determine the potentials <fi and A z as a function of x and y under a certain 
boundary condition, we obtain the distributions g^o for each species, and vice versa. 

Here we note restriction of the method described above. In general, a stationary equi- 
librium can have various velocity distributions as long as the pressure balance is achieved. 
However, Equations ^ look as if there is a one-to-one correspondence between the field con- 
figuration and the distributions gj$. This disagreement comes from the assumption that the 
deviation of the velocity distribution from the Maxwellian can be expanded by the Hermite 
polynomials. Therefore, among a lot of possible equilibria, the stationary solutions derived 
above cover a part of them. 



III. MAXWELLIAN EQUILIBRIUM 



In this section, we consider applications of the stationary solutions derived above. Most 
of familiar stationary solutions to the Vlasov-Maxwell system are based on the classical 
Maxwell-Boltzmann distribution. For example, the Harris sheet equilibrium [18] and the 
Bennet pinch equilibrium [19j are well-known one-dimensional equilibria. Our solutions 



include these equilibria. We reproduce the Harris sheet equilibrium in |III A| and the Bennett 
pinch equilibrium in III B| by using our method. 
At first, we assume that gj t0 is expressed as 



9j,o = exp 



rrijVjC 



(10) 



where Uj is a constant with the dimension of velocity. Then, using relation (JTl), we obtain 



9j, 



u , 



ni \ v 



exp 



2q j u j A z 
nijv'jC 



(11) 



Substitution of this expression into Equation (J5| yields 



fi 



oo 1 

Y- 



n=0 v 



U, 



H„ (v z /vj) 



n.; 



71-3/2^3 



exp 



V 2 + V 2 + V 2 

x y z 



rrijV 2 



UjA z 



(12) 



By performing the summation with respect to n in this equation through the relation (A3) 
we obtain a shifted Maxwell-Boltzmann distribution; 

\2 



fi 



n 



7j-3/2-?.3 
3 



exp 



4 + v 2 + (v z 



V) 



2qj_ 
rrijV 2 



UjA, 



(13) 



From this expression, it is found that Uj is interpreted as the bulk velocity of species j. 
On the other hand, substitution of Equation ( 10 ) into Equation ^ leads to 



d 2 



+ 



d 2 



dx 2 dy 



d 2 A 7 d 2 A 2 
+ 









jn^iexp 


f- 2? ' 2 ( 






rriivf \ 





(14a) 



+n e q e exp 



2q e 



U e ,A, 



dx 2 



dy 2 



47T 

c 



jn^iMi exp 


2qi ( UiA z \ 
rriiV 2 V c J 


+n e q e u e exp 


m e vi \ 



(14b) 



Up.A, 



In order to obtain the electromagnetic fields as a function of x and y, we need to solve 



these equations under a certain boundary condition. Although Equations (14) seem to be 



complicated, we can derive some analytic solutions under appropriate assumptions. 
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In the following derivation, we assume that the plasma comprises of electrons and ions 
with same charge but the opposite sign (g^ = — q e = e) and satisfies the charge neutrality 
(n, = n e = n ) and that the electric field strength is sufficiently small (0 = 0). Then, 



Equation ( 14b ) leads to 



U.; 



rriivf 



m e v, 



(15) 



which implies that electrons and ions move in opposite directions. 



A. Harris sheet equilibrium 



The Harris sheet equilibrium [18J is a well-known Vlasov-Maxwell equilibrium, which is 
widely used in various analytical and numerical studies of collisionless plasmas. For example, 
this equilibrium provides the initial setup for investigations of the collisionless tearing mode 
instabilities. 

In order to derive the Harris sheet, we assume that the plasma is homogeneous in y- 



direction. Combining with the relation (15), Equation (14c) reduces to the following second- 
order ordinary differential equation; 

d 2 A z 4Trn e(ui-u e ) (2euiA z \ 

TT = ex P 2~ ■ ( 16 ) 

dx z c \ rriivfc J 

This equation has a solution in the form of 

A z = —B L\n[cosh(x/L)}, (17) 
where B and L are constants satisfying the following relations; 

B 2 = 47m (miV* + m e v 2 e ), (18) 

and 

B L = ^ = (19) 

eui eu P 



The expression (17) gives a sheared magnetic field; 

B y = B Q tanh(x/L), (20) 
and a distribution function in the form of 

A K (21) 

From these expressions, Bq is interpreted as the asymptotic field in the limit of x — > oo and 
L the scale-length. 



B. Bennett pinch equilibrium 



The Bennett pinch [19J is a Vlasov-Maxwell equilibrium in the cylindrical configuration. 



In cylindrical symmetry, Equation ( 14c ) reduces to the following second-order ordinary dif- 
ferential equation. 



1 d ( dA 



47m e(ui — u e ) ( ItUiA 
exp 



rriivfc 



r dr \ dr 

where r = 'x 2 + y 2 . This equation has a solution in the form of 

A z = -B L\n [l+r 2 /(4L 



(22) 



(23) 



Also in this case, two constants satisfy the relations (18) and (19). Then we obtain a 
self-consistent cylindrical configuration of the plasma; 

1 



fj oc 



[l + r 2 /(4L 2 )] 2 ' 
and it generates a magnetic field surrounding itself; 



(24) 



B, 



2B Q Ly 2B Lx 

A; 



(25) 



r 2 + 4L 2 ' y r 2 + AL 2 ' 

In this case, Bq and L are interpreted as the scale of the magnetic field strength and the 
scale of the current tube, respectively. 



IV. NON-MAXWELLIAN EQUILIBRIUM 

In the previous section, we have dealt with one- dimensional Vlasov-Maxwell equilibria 
based on a shifted Maxwell-Boltzmann distribution. Furthermore, by applying our method, 
we can construct two-dimensional equilibria with a distribution function deviating from the 
Maxwell-Boltzmann distribution as follows. 

First, we assume the charge neutrality and a small electric field strength again. Then, 



Equation (9a) implies that ions and electrons have the same distribution; 

9i,o = 9e,o- (26) 
Furthermore, we assume that A z takes the following form, 

A z = -B L[cos(x/L) + cos(y/L)], (27) 



which provides a sinusoidal magnetic field as 



B x = B sm(y/ L), B y = -B sin(x/ L) . 



(28) 



Substituting Equations (26) and (27) into Equation (9b), we obtain the following first-order 



ordinary differential equation; 



dgt. 



i.O 



dA z 27in L 2 (mivf + m e v 2 ) 
Integration of this equation gives 

A 2 + C 

9i,0 — 9e,0 



(29) 



AimnL 2 (rriivf + m e v 2 



(30) 



where C is a constant of integration. Substitution of this expression into Equation ([7]) leads 
to 



9j,i 



9j,2 



rrijVjcA z 



AiTq j noL 2 (mivf + m e v 2 ) 
m 2 VjC 2 



167Tq 2 n L 2 (miV 2 + m e v 2 ) 
g j>n = for n > 3. 



(31) 

(32) 
(33) 



As a result, we obtain the following distribution function after simple algebraic manipula- 
tions; 

2 



fi 



1 



47T 5//2 L 2 (mjf 2 + m e v 2 )v 3 



A z + 



rrijC 



+ A 2 



exp 



v 2 + v 2 + v 2 

x 1 y 1 z 



(34) 



Here, the explicit expressions for H , Hi, and H 2 (A4) have been used and A is a constant 
defined as 

2 7rijV 2 c 

An — C ~ f; , (35) 



1 



in order that the distribution function (34) takes positive values at arbitrary points in the 



phase space. This distribution satisfies the Vlasov-Maxwell system ([T]) - Q exactly. 

Figure [T] illustrates the thus derived configuration of equilibrium. The gray scale and 
the arrows represent the density distribution gjp of electrons or ions in arbitrary units and 
the magnetic field, respectively. We can see that the current filaments lie along the z-axis 
and generate the magnetic fields around themselves. Each filament is surrounded by four 
filaments that carry anti-parallel currents. 



Here, we focus on the velocity distributions at two points. One is the center of a filament 
(refer to as the O-point) and the other is a point between the filaments (refer to as the 
X-point). The magnetic field at the X-point is sheared one like the Harris sheet equilibrium. 
Figure [2] shows the velocity distributions of electrons at the X and O-point, respectively. 
The values of parameters are as follows; Bq = 0.04m e u> e /e, where u e is the electron plasma 
frequency, A = 0, and v e = 0.1c. At the X-point, we see the double-peak distribution. 
On the other hand, a distribution similar to the shifted Maxwellian one is achieved at the 
O-point as long as eB L > m e cv e holds. 

Finally, we mention the velocity distribution of ions. Because of the large mass 



(mi ~ 2000m e ) and the small thermal velocity (vi <C v e ) according to Equation (35), the 



distribution (34) has a sharp double-peak like two delta functions and becomes almost in- 
dependent of x and y. 



V. CONCLUSIONS 



In this paper, we have developed a novel method to construct stationary solutions of the 
Vlasov-Maxwell system. This investigation must lead to great insight into the properties 
of self-consistent plasma configurations. By applying our method, some of previous studies 
are reproduced and a new two-dimensional equilibrium is proposed. In particular, it is 
found that a sheared magnetic field like that in the Harris sheet equilibrium is generated. 
Although the stability of the equilibrium configuration needs to be investigated, it may 
provide an initial setup for investigations of the collisionless magnetic reconnection that 
takes into account three-dimensional effects. 
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APPENDIX A: PROPERTIES OF HERMITE POLYNOMIALS 



In this section, we review the properties of the Hermite polynomials used in this paper 
(see, for example, [22] )• The nth-order Hermite polynomial H n (x) is defined as the solutions 
of the Hermite's differential equation; 

-ry - 2x— - + 2nH n = 0, Al 
ax z ax 

The set H n (x) forms an orthogonal basis in the function space. In other words, the Hermite 
polynomials satisfy the orthogonality relation; 

/oo 
H m (x)H n (x)e~ x2 dx = 2 n yfrn\5 m , n , (A2) 
■oo 

where 5 m ^ n represents Kronecker's delta. The Hermite polynomials are also given by the 
generating function as 

e -t* + 2 Xt = * Hn (x). (A3) 

n=0 

From this equation, we can obtain the expression for the Hermite polynomials. The first 
three are 

H = 1, H 1 = 2x, H 2 = 4x 2 - 2. (A4) 
The derivative of H n is expressed in terms of H n -\ as 

^ = 2nH n ^. (A5) 
ax 
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